The dynamically asymmetric SQUID: Miinchhausen effect 
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We report on a complex zero-temperature decay channel of a classical object in a metastable state 
coupled to a quantum degree of freedom. This setting can be realized in a dc-SQUID where both 
Josephson-j unctions have identical critical currents Ic but feature strongly asymmetric dynamical 
parameters; more precisely, selecting both parameters C and 1/R adequately large for one and small 
for the other junction makes the first junction behave essentially classical but lets quantum effects 
be present for the second one. The decay process is initiated by the tunneling of the quantum 
junction, which distorts the trapping potential of the classical junction; the metastable state of the 
latter then becomes unstable if the distortion is large enough. We present the dynamical phase 
diagram of this system providing the dependence of this decay channel on the external bias current 
/ and on the coupling strength between the two junctions, determined by the loop inductance L. 

PACS numbers: 85.25.Dq, 74.50.+r 
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I. INTRODUCTION 

It is a basic feature of the classical world that a mas- 
sive object residing in a metastable potential weU cannot 
decay and, at zero temperature, is determined to reside 
in the well. There are, however, situations, where this 
doctrine is no longer adequate, namely if the system con- 
sists of two parts, a classical and a smaller quantum one. 
In that decay process is possible under certain 

circumstances: The decay sequence is initiated by the 
tunneling of the quantum degree of freedom; through its 
coupling to the classical part the trapping potential is 
distorted, eventually turning the metastable state into 
an unstable one, provided that the distortion is large 
enough. We name this decay process the 'Miinchhausen 
effect' after the famous baron telling the story of pulling 
himself out of a swamp by his own hair. 

The above situation can be implemented experimen- 
tally in a dynamically asymmetric dc-SQUID; it is this 
specific realization we will study in detail in this pa- 
per, see Fig. [T] The dynamical degrees of freedom in 
the dc-SQUID are the gauge-invariant phase differences 
(pj, j = 1, 2, across the two Josephson junctions. The po- 
tential energy (of a single Josephson junction) is given by 
Vj = Ej(l—cosipj), j — 1,2, involving the Josephson en- 
ergy Ej = $o^c/27rc (with the flux unit $0 = hc/2e and 
the critical current Ic of the junction). The kinetic en- 
ergy reads Tj = (fi/2e)^Cj(^|/2, where the capacitances 
Cj assume the role of effective masses. Hence, a SQUID 
featuring two Josephson junctions with equal critical cur- 
rent Ic but adequately chosen and strongly asymmetric 
capacitances, one big and one small (Ci 3> C2), effec- 
tively provides us with a classical and a quantum degree 
of freedom. The fabrication of large, classical junctions 
is easily achieved today; however, this is not the case for 
small quantum junctions whose realization is more dif- 
ficult. Nevertheless, small junctions exhibiting quantum 
tunneling^'^ and quantum coherence^i^ can be routinely 



"classical" 




FIG. 1: Schematics of the dynamically asymmetric dc- 
SQUID. Two Josephson junctions with equal critical current 
Ic but strongly asymmetric capacitances (Ci ^ C2) and re- 
sistances Rj are integrated in a biased (current 7) supercon- 
ducting loop with symmetric inductance L. 



fabricated today. 

The decay process of the biased dynamically asymmet- 
ric SQUID proceeds in the following manner: The bias /, 
leading to a term oc —I{(pi +'P2) in the potential, turns a 
stable state of the washboard potential into a metastable 
one. As junction 1 features a large capacitance, we as- 
sume its dynamics to be strictly classical. If the bias / is 
large enough, a (imaginary time) decay process involving 
only the quantum junction (at constant tpi) is enabled. 
This phase slip leads to the entry of magnetic flux into 
the ring. Given the inductive coupling oc 1/L (inductance 
L), the current through the classical junction is enhanced 
and it may eventually become overcritical, thus decaying 
via a classical real time trajectory. 

In the following, we specify the setup (Sec. [n|. In 
Sec. mil we describe the decay process in more detail 
and present the resulting dynamical phase diagram. We 
finish with a few concluding remarks in Sec. IIVI 
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II. SETUP 

We start from the capacitively shunted junction model 
(CSJ), where the dc-SQUID, biased with a current /, is 
described by the Lagrangian 



quantum junction are incoherent^ii. We will analyze 
this situation in detail in the next section; the obtained 
results are also relevant for other choices of parameters, 
c.f. below. 
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here, we have assumed that the inductance L of the 
SQUID is symmetrically distributed. The Lagrangian 
C generates the equations of motion 



(2) 



with the 'masses' 



^(jCj /2ttcIc and where we have 



added the dissipative terms r/jipj with the damping pa- 
rameters rjj — $o/27rc/c-Rj oc i/Rj, the normal resis- 
tances of the junctions; the potential (illustrated in Fig. 
[2]) is given by 



C0Siy9i — C0S(p2 
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with the dimensionless current i — I /2Ic, the coupling 
constant k = ^o/2t:cIcL and where energy is measured 
in units of Ej. 

Residing in a symmetric {(pi — ^2 — arcsini) 
metastable state of the potential v{(fi,ip2) at finite bias 
current i, the classical version of the system described 
through Eq. ([2]) cannot decay at zero temperature. Here, 
we are investigating the case where junction one, featur- 
ing a large capacitance Ci , is assumed to behave strictly 
classical, whereas the dynamics of junction 2 is char- 
acterized by large quantum fluctuations.^ This can be 
achieved through a suitable choice of parameters, i.e., 
small capacitance C2, while keeping Ej > Ec2 — e^/2C2 
such that we remain in a quasi-classical regime. 

In the system under consideration, different scenarios 
can arise depending on the strength of the dissipation as 
quantified by the dimensionless damping parameter 



(4) 



where, at z = 0, the plasma frequency hujpj — 
{SEjEcjY^^- The simplest case is the overdamped situa- 
tion ai,a2 > 1, where the dynamics of the classical junc- 
tion is viscous and both relaxation and tunneling of the 



^ Hence, a description via an equation of motion, Eq. ((2]l, is not 
suitable for the quantum junction and dissipative effects have to 
be included through a path integral formalism. 



III. DECAY PROCESS AND PHASE DIAGRAM 

In the following, we determine for which currents i and 
coupling constants k a zero temperature decay of a sym- 
metric metastable state {ipi = ip2 = arcsini, up to an 
arbitrary multiple of 27r) is allowed in the interferometer 
potential w(<pi, <^2)- The result is displayed in a dynami- 
cal phase diagram in the j-fc-plane, see Fig. [31 where the 
critical line ic(fc) separates regions where this decay is 
prohibited (localized state) from regions where it is al- 
lowed (delocalized). 

Assuming junction 1 to behave strictly classical, (i.e. 
considering the limit of very large Ci), a quantum decay 
of the SQUID in a metastable state can only occur at 
fixed i.e. through an imaginary-time trajectory of ip2 
in the effective potential Wcff(V2) — v{^i = const., (^2)- 
We will adopt this approximation for all tunneling pro- 
cesses throughout the discussion. 

In the case ai, a2 > the preparation in the assumed 
symmetric initial state is straightforward: A SQUID is 
cooled close to zero temperature at zero external mag- 
netic field before a bias current is ramped to the desired 
value i. The ramp needs to be fast enough, such that the 
final current is reached before the decay of the quantum 
junction; any initial potential energy is dissipated as the 
phases relax to the bottom of the metastable well. 

For undercritical currents i < 1 the described ini- 
tial state, fi — f2 = arcsini, is stable against a de- 
cay involving the classical junction; however it is, for 
k < i/(7r — arcsini), unstable with respect to a quantum 
decay of ip2 since the minimum of Ves{^2) near ip2 ~ 2tt 
is lowered below the initial one and one or more phase 
slips of ip2 are possible. In order to determine whether 
the SQUID will remain in a stable or enter a finite volt- 
age state at given i, k we proceed in two steps: First, we 
have to determine where the successive quantum tunnel- 
ing of ip2 comes to a halt, i.e. which side minimum is 
quantum-stable. This is equivalent to finding the global 
minimum of the effective potential Vcsi^2)- Fixing the 
phase across the classical junction ipi = arcsinz, one im- 
mediately sees that the quantum-stable minimum is the 
one near ip2 ~ 27rn if n is the largest integer such that 



k < 



(2n — l)7r — arcsini 



(5) 



The sequence of phase slips of the quantum junction 
leads to an accumulation of magnetic flux in the SQUID 
loop, inducing a screening current. Consequently, the 
current through junction 2 is reduced whilst that through 
junction 1 is increased. The magnitude of the induced 
current depends strongly on the coupling constant fc; for 
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given n, it is only large enough to drive the classical junc- 
tion overcritical if 



1 



(2n - 1/2)tt + arcsin(2i - 1) 



(6) 



If this is the case, the SQUID enters a finite voltage state 
where quantum tunneling of junction 2 (an approximate 
flux unit enters the loop) and classical relaxation of junc- 
tion 1 (flux leaves the loop) interchange sequentially (see 
Fig. [2]). If /c < fcc,n(«), the SQUID resides in a localized 
state and a further increase in i is necessary to drive the 
system unstable. The two conditions Eqs. (O and ([6]) 




FIG. 2: Illustration of the decay sequence of the dynami- 
cally asymmetric SQUID [i = 0.5, k = 0.04). The initial 
metastable well is unstable w.r.t. the macroscopic quantum 
tunneling of the small junction 2. A continuous sequence of 
phase slips takes the system to a state which is classically un- 
stable. In the following, the (classical) relaxation of the large 
junction 1 and the quantum decay of junction 2 alternate and 
lead to a finite voltage state of the SQUID. 

generate a web of crossing lines in the i, fc-plane, deter- 
mining the critical line ic(fc) marking the dynamic tran- 
sition from a localized to a delocalized state (Fig. [3]). In 
the limit fc — > 0, where the cosine in the potential Eq. ([3]) 
becomes a small correction to the parabola, the critical 
line ic{k) approaches 1/2; this indicates that all current 
is redirected through junction 1 and dclocalization takes 
place at / = Ic, the critical current of a single junction. 
The simple arguments above have to be refined in 
order to obtain the precise location of the critical line 
ic{k). First, the condition of classical stability is but 
the standard determination of the critical current of a 
dc-SQUID's asymmetric minimuni^. In our case, where 
classical stability along the (/92-direction is guaranteed, 
the relevant set of equations is given by 

sin((^5') = i-kc,n{'fi ~'f2)> (7) 
Sm{<f'^) = i + kc,n{<Pi -'f2)^ (8) 
cos COS (^2 — —fcc.n (cos + cos 1^2); (9) 

and has to be solved (numerically) for kc.n{i) and 
(^"21 the coordinates of the true minima near ipi = 



arcsini,<y92 = 27rn. Eq. ([6]) is an approximate solution 
to Eqs. Q-® in the hmit of fc ^ 1. Second, as tunnel- 
ing of the quantum junction might be enabled only after 
the relaxation of the classical junction to a minimum, 
condition Eq. (O has to be corrected to 



kc,n{i) 



{In - l)7r-(^5'- 



(10) 



taking into account the change of the effective potential 
Va?i{}P'i) upon a change in ip\. The exact numerical so- 
lutions of Eqs. Q-lIll) and ([T0|) are shown in the inset 
of Fig. [31 where the approximate solution is seen to be 
rather precise. 
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FIG. 3: Dynamical phase diagram of the dynamically asym- 
metric dc-SQUID. The critical line (solid) separates regions 
in the i-fc-plane where the SQUID enters a continuous finite 
voltage state (delocalized) from regions where the phases re- 
main localized. On the lower-bias side of the critical line, 
the dotted lines mark the entry of a additional flux into the 
SQUID loop. On the higher-bias side of the critical line, the 
dashed lines show how many flux units are at least needed to 
render the SQUID unstable. The inset shows a close up on 
the critical line and additionally displays the exact numerical 
solution (dashed). 

A few remarks on the above results are in place. First, 
the setup as described above may not be suitable for 
experimental investigations as an overdamped quantum 
junction suffers from strongly suppressed tunneling rates. 
However, the analogous results can be obtained by us- 
ing two underdamped junctions and ramping the cur- 
rent sufficiently slowly, such that the junctions have dis- 
sipated all energy before crossing any line (solid or dot- 
ted) in the phase diagram Fig. [S] The hallmark of the 
"Miinchhausen-effect" in this case is the two different 
types of decay that finally initiate a running state: If 
the critical line is crossed on a part with positive slope, 
the quantum stable minimum is classically unstable im- 
mediately and a finite voltage state is initiated by the 
last quantum tunneling process. On the other hand, if 
the critical line is crossed on a part with negative slope, 
the quantum-stable minimum is still classically stable af- 
ter the tunneling and needs to be turned overcritical by 
increasing the bias. Then, the finite voltage state is initi- 
ated by a purely classical decay. The difference between 
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the two types of decay should be visible upon analyzing 
the decay histograms of multiple measurements, being 
broad for a quantum decay but narrow for a classical 
process. 

Second, one has to take into account that for realistic 
systems, where the capacitance asymmetry is introduced 
by a capacitive shunt^, the maximum asymmetry is lim- 
ited. Corrections due to the finite mass of the classical 
junction might appear in the form of two-dimensional 
quantum tunneling, where junction 1 takes part in the 
tunneling process^^. 



tions behaves fully classical, provided the second junction 
shows quantum behavior. The "Miinchhausen-decay" in- 
volves tunneling of the quantum junction, where mag- 
netic flux accumulates in the superconducting ring and 
eventually redirects enough current through the classical 
junction as to drive it overcritical. We thank A. Ustinov 
and A. Wallraff for fruitful discussions. 



IV. CONCLUSIONS 



We have shown that a dc- SQUID can decay out of 
a symmetric metastable state, even if one of the junc- 
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